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ABSTRACT 

Context. The formation of massive stars goes through phases which remain heavily obscured until the object is well on the main se- 
quence. The identification of massive YSOs in different evolutionary phases is therefore particularly difficult, and requires a statistical 
approach with large samples of candidate objects to determine the observational signatures of these different phases. 
Aims. A mandatory activity in this context is the identification and characterisation of all the phases that a massive forming YSO 
undergoes. It is of particular interest to verify the observability of the phase in which the object is rapidly accreting while not yet 
igniting the fusion of hydrogen that marks the arrival on the ZAMS. 

Methods. One of the candidate prototypical objects for this phase is Moll60/IRAS23385+6053, which has been the subject of de- 
tailed studies which confirmed and strengthened the possibility that this massive YSO may be in a pre-Hot Core stage. We further 
investigate this issue by means of Spitzer imaging and spectroscopy in the 5-70 yum range. 

Results. The dense core of Moll60/IRAS23385+6053, which up to now had only been detected at submillimeter and millimeter 
wavelenghts, with only upper limits below 20 //m, has been revealed for the first time at 24 and 70 fim by Spitzer. These observations 
confirm the earlier assumptions that this object is dominant at far-IR wavelengths. The complete 24 /jm-3.4 mm continuum cannot 
be fitted with a standard model of a Zero-Age Main-Sequence (ZAMS) star embedded in an envelope. A simple greybody fit yields a 
mass of 220 Mq. The luminosity is slightly in excess of 3000 Lg, which is a factor of 5 less than previous estimates when only IRAS 
fluxes were available between 20 and 100 yum. The source is under-luminous by the same factor with respect to UCHii regions or Hot- 
Cores of similar circumstellar mass, and simple models show that this is compatible with an earlier evolutionary stage. Spectroscopy 
between 5-40 pm shows that the physical conditions are typical of a photo-dissociated or photo-ionised region, where the required 
UV illumination may be provided by some of the other sources revealed at A <24pm in the same star forming region, and which can 
be plausibly modeled as moderately embedded intermediate-mass ZAMS stars. 

Conclusions. Our results strengthen the suggestion that the central core in Moll60/IRAS23385+6053 is a massive YSO actively 
accreting from its circumstellar envelope and which did not yet begin hydrogen fusion. 
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1. Introduction 

The identification and the characterisation of all the different 
phases that an intermediate and high mass forming object under- 
goes during its approach toward the Main Sequence (MS) has 
received increasing attention in recent years. Besides its impor- 
tance to complete our understanding of the star formation pro- 
cess in general, this research has interesting consequences in a 
broader context. For example, the number of OB ZAMS stars 
indirectly estimated using a variety of tracers like the radio con- 
tinuum or the Hff emission, is widely used to estimate the star 
formation rate and -efficiency, critical parameters to measure the 
history of star formation in our Galaxy as well as in external 
galaxies (McKee & Williams [19971 Kennicutt |T998| . 

The identification of a "protostar" for intermediate and high- 
mass objects is more difficult than for those of low-mass, be- 



cause massive objects reach the ZAMS while still accreting ma- 
terial. Up to that moment they remain deeply embedded in their 
natal molecular cloud, and unlike their lower-mass counterparts 
they lack a visible pre-MS phase. With accretion rates lower than 
10"^ Mq yr"' the existence of this classical pre-MS phase is lim- 
ited to objects of mass M<8 Mq. Accretion, however, plays a 
very important role raising this mass limit and allowing the ar- 
rival on the ZAMS at higher stellar masses if accretion proceeds 
at much higher rates; the Young Stellar Object (YSO) joins the 
ZAMS at M= 16 Mq for M= lO""* Mq yi--' (Palla & Stahler l992] l, 
and well beyond 20 M© for even higher rates. Such rates seem in- 
deed plausible, given the mass loss rates measured from outflows 
emanating from massive YSOs (Zhang et al. |200U Beuther et al. 
2002), and are also predicted by theoretical models (McKee & 
Tan 2003) where accretion rates increase with the mass of the 
forming core. 
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We recently completed a preliminary study of the Spectral 
Energy Distributions (SEDs) of a sample of massive YSOs in 
various evolutionary stages (Molinari et al. 2008) which sug- 
gests that objects more massive than 8 may indeed be 
observable in a pre-ZAMS phase as pure far-infrared (unde- 
tected below 20 fim) objects underluminous with respect to Hot- 
Cores or Ultra Compact (UC) Hii regions for a given enve- 
lope mass. Moll60/IRAS23385+6053 clearly manifests its dif- 
ference in multiwavelength appearance with respect to known 
Hot-Cores or UCHii regions, and is perhaps a template for 
this class of massive analogues of classical Class objects 
(Molinari et al. 1998b). The strong millimetric peak, which 
appears to harbour the driving source of a molecular outflow, 
is anti-correlated with the complex emission morphology vis- 
ible at mid-IR wavelengths (Molinari et al. 1998b) and with 
the extended radio continuum emission found in the region 
(Molinari et al. 120021 . Detailed spectroscopic investigation con- 
firmed that from a chemistry viewpoint, the millimeter core of 
Mo1160/IRAS23385-h6053 appears to be on t he verg e of turn- 
ing into a Hot-Core (Thompson & Macdonald 2003") and, with 
T<50K (Fontani et al. 2004), is considerably colder than typical 
Hot-Cores (HC) of similar envelope mass. 

The suggestion that bolometric luminosity and SED mor- 
phology may be promising evolutionary indicators is based on 
the coupled analysis of mid-IR and submillimetric data. For the 
study of large samples in these wavelength ranges, where sur- 
vey data are available, the accessible spatial resolutions of 10- 
20" are not sufficient to pinpoint each and every emitter in the 
clustered environments where massive YSOs are found, but are 
enough to recognize systematic differences between classes of 
objects. Similar resolutions in the far-infrared, where they are 
most crucial since it is where the SEDs peak, will not be ac- 
cessible until the Herschel satellite is available; in the mean- 
time, an important confirmation of the scenario proposed by 
Molinari et al. (2008) may be provided by the Spitzer satel- 
lite. Stecklum et al. (.2005) presented MIPS continuum data of 
Moll60/IRAS23385-i-6053 and concluded that this source was 
not a high-mass protostar. We here reconsider the same MIPS 
data and provide an alternative and more plausible intepretation 
which is more consistent with the independent evidence of the 
pre-HC nature of this object and the more general picture of the 
luminosity evolution of massive YSOs. Additional Spitzer IRS 
spectroscopy is also presented, and the implications for the star 
formation activity of the entire region are discussed. 



2. Observations and Data Analysis 

Data for source Moll60/IRAS23385-i-6053 were acquired with 
the Spitzer satellite in two Guaranteed Time programs. The first 
program used the MIPS instrument (Rieke et al. 2004^ to acquire 
photometric images at 24 and 70 //m. The 128 x 128 Si:As array 
at 24 //m has a pixel size of 2." 49 x 2."60 and a field-of-view 
(FOV) of 5. '4 X 5. '4, while the 70 pm amy consists of 32 x 32 
Ge:Ga matrix with a pixel size of 5."24 x 5"33 and a FOV of 2. 
'7 X 1.'35. The data were taken with an exposure time of 3 sec, 
and were repeated twice at 24 pva and 3 times at 70 /im. One 
cycle of the basic 24 pm photometry observation consisted of 
14 offset images of the field combined with classical techniques 
of chopping (along-scan mirror motion) and dithering. One cy- 
cle of 70 jum observations produces, with the same method, 12 
images of the field. Also an estimate of background-foreground 
contribution is made acquiring an image off-source with chop- 
ping technique. 



The second program contains spectroscopy mapping with 
the IRS instrument (Houck et al. 2004) using three of the avail- 
able modules. The "Long High" (LH) module covers a range 
between 19 and 37 pm. at a spectral resolution of R=600 with a 
FOV of 22."3 X 1 1 ." 1 and a pixel size of 4."5. The "Short High" 
(SH) module covers with similar spectral resolution the range 
between 10 and 19.5 pm with a FOV of 1 1."3 x 4."7 and a pixel 
size of 2."3. The "Short Low" (SL) slits were used to cover the 
range between 5.6 and 14.5 pm. The SL module is composed of 
two sub-modules at different wavelengths; one covers the range 
between 7.4 and 14.5 pm at R=16, with a FOV of 57"x3."7 and a 
pixel size of l."8, and the second covers at R=8 a range between 
5.2 and 8.7 pm, with a FOV of 57"x3."6 and a pixel size of 1. 
"8. Observations consisted of one integration with a ramp dura- 
tion of 6 seconds. The area around the position of the millimeter 
emission peak was covered with a small 3x2 map. 

For both programs we downloaded the Basic Calibrated Data 
products from the Spitzer archive. Pipeline versions used for au- 
tomatic processing were S.l 1.0.2 for MIPS data, and S 15.3.0 for 
IRS. 

2.1. Source photometry: Near-IR to millimeter 

Fig. [Bi shows the Mo1160/IRAS23385h-6053 field at 24 pm 
(greyscale image) and at 70 pm (contours). The extraction and 
photometry of point sources from the 24 pm image was at- 
tempted using several publicly available packages for automatic 
source detection and photometry (SExtractor, IRAF, MOPEX). 
As is apparent from the image, the various peaks that are visible 
sit on a complex patch of extended emission which is strongly 
variable on small and large scales. This poses severe problems 
for reliable source detection and photometry since the various 
packages we tried assume either a constant background, or esti- 
mate one by using a box of fixed size across the image. Sources 
like A in the figure, for example, are undetected with SExtractor 
and also with MOPEX if the box for estimating the background 
is too large. Using a small box would solve this, but on the 
other hand leads one to significantly understimate the fluxes for 
sources like N or E which instead lie on a less variable and 
broader patch of extended emission; discrepancies in the fluxes 
of up to a factor of 5 are found for some objects using different 
packages. For this reason we decided to visually identify com- 
pact peaks in the map, and to manually perform the photometry 
at 24 pm with custom procedures using the GILDAS package, 
optimising by hand the area for the background estimate. 

At 70 pm the situation is less complex because due to 
the lower spatial resolution the various sources in the area are 
blended in a large patch of emission; among the sources detected 
at 24 pm source A is the one that more clearly stands out above 
this plateau (TablelTJ. A slight offset is apparent between the po- 
sition of the 24 and IQpm peaks; we will comment on this below. 
The size of this peak was estimated analyzing cross-cut bright- 
ness profiles; this information was transferred to the AIPS task 
JMFIT to fit a composite 2D-Gaussian and a plateau resulting in 
an integrated flux density of 65 Jy for the compact central peak 
coincident with source A. The integrated intensity of the entire 
70 pm emission structure, after subtraction of a constant emis- 
sion level clearly detectable above the noise at the borders of the 
image, is ~ 350 400 Jy, compatible with the 60 pm IRAS flux 
density. The 70 pm flux density that can be assigned to the cen- 
tral massive core coincident with source A is then much lower 
than the 60 pm IRAS flux density, contrary to the assumptions 
made by Fontani et al. (2004). In that work we carried out a far- 
IR extrapolation of the mid-IR emission of the bright patches 
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Fig. 1. a) (Left) MIPS 24 yum grey-scale image of the Moll60/IRAS23385+6053 field, with superimposed the contours of the 70 /urn 
emission; the letters mark the identified sources, b) (Right) 2.2 fj.m image of the same region, at approximately the same scale; solid 
countours represent the 24 fim continuum emission, while the dashed contours show the 3.6 cm thermal free-free emission. 



Table 1. Photometry of sources in the field 



Source 


Q'(J2000) 


<5(J2000) 


Fb 






F« 




F" 

6.75^/(1 


I5fjm 




FvOf/zii 


FgSOfim 










mJy 


mJy 


mJy 


mJy 


mJy 


Jy 


Jy 


Jy 


Jy 


Jy 


A 


23:40:54.65 


+61 


10:27.6 
















0.14 


55 


2.2 


B 


23:40:50.42 


+61 


10:37.4 






0.8 


1.5 


1.7 


0.2 


0.14 


0.16 


<25 


<0.3 


C 


23:40:50.25 


+61 


10:27.6 






1.5 


6.1 


10.6 


<0.02 


0.6 


0.6 


<25 


<0.3 


D 


23:40:57.11 


+61 


10:41.0 






0.07 


0.35 


0.7 


0.2 


0.1 


0.3 


<25 


<0.4 


E 


23:40:56.35 


+61 


10:30.0 






0.07 


0.16 


0.6 


0.3 


0.2 


1.0 


<35 


<0.8 


F 


23:40:55.50 


+61 


10:49.0 












0.26 


0.2 


0.5 


12 


<0.5 


G 


23:40:52.62 


+61 


11:09.2 






0.1 


0.06 


0.2 


e 




0.13 


<30 


<0.5 




fl 23:40:58.97 


+61 


10:27.9 


0.1 


04 


0.07 


0.5 


0.5 












H* 


•! 2 23:40:59.04 


+61 


10:33.3 


0.05 


0.2 


0.07 


0.1 


0.1 


<0.01 


<0.15 


0.14 


<20 


<0.2 




^3 23:40:59.11 


+61 


10:43.8 






0.8 


0.8 


0.6 












I 


23:40:57.62 


+61 


11:08.0 


0.2 


0.7 


6.4 


6.6 


5.1 


0.24 


0.1 


0.08 


c 


<0.2 


L 


23:40:49.23 


+61 


11:23.9 






0.64 


2.0 


5.5 


e 


e 


0.3 


c 


<0.1 


M 


23:40:47.45 


+61 


10:57.6 












0.04 


0.01 


0.33 


<15 


<0.2 


N 


23:40:41.44 


+61 


10:28.1 






0.5 


1.5 


2.4 


e 


e 


0.4 




_d 


P 


23:40:49.23 


+61 


10:11.6 






0.4 


1.7 


4.9 


0.02 


0.03 


0.1 


<15 


<0.4 


R 


23:40:48.55 


+61 


10:15.9 






1.0 


1.3 


2.1 


f 


e 


0.08 


<15 


<0.4 



° from ISOCAM images (Molinari et al. il998bi ) 

* Although the positional association of other K, sources to the 24 yum peaks is very good (see Fig.lJJ)), it seems likely that source H, since it is 
much shallower than the other 24 fim peaks, could be associated to one of the three nearest K,, sources visible in its close proximity. 

This source is not covered by the 70 fim image. 

This source is not covered by the 850 /jm image. 
*■ This source is not covered by the ISOCAM images. 



revealed in the ISOCAM images (Mohnari et al. |1998bl ), which 
are also visible in the present 24 /vm image, and concluded that 
most of the IRAS 60 and 100 fim flux densities could be assigned 
to the central core. The MIPS observations, however, provide a 
direct measurement and show that only a minor fraction of the 
total A <10 fim flux density in the region can be assigned to the 
central core. 

In addition to source A, which is clearly dominant at 70 fim, 
we detect a fainter peak which barely emerges from the plateau 
in correspondence to source F; as for source A, the integrated 
flux is estimated using the AIPS JMFIT task fitting a peak- 



plateau combination. Another secondary peak is visible both at 
24 and 70 yum about 35" South/South-East of A, but is consid- 
erably shallower and broader than the other sources and is not 
considered in this analysis. As for the other sources detected at 
24 jum, we assigned the local values of the 70 fim emission as 
upper limits. 

The 24-70 jum spectral region alone, although very important 
as the SED rises toward the peak for typical YSOs envelopes, 
is not sufficient for a reliable comparison with SED models. 
Leaving aside IRAS whose beam at 100 yum encompasses more 
than the entire area shown in Fig.[T^, the longer wavelength in- 
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Fig. 2. Full IRS spectrum (top panel) extracted using CUBISM at the position of the central millimeter core (source A). The bottom 
panel is an enlargement where this time the SH module spectrum is resampled at the resolution of the SL module; its rescaling 
by a factor of 4.2 brings it to an almost perfect overlap with the SL spectrum, and joins very well with the LH spectrum at longer 
wavelength. 



formation is extracted from submillimeter and millimeter im- 
ages (Molinari et al. 120021 1. Similarly to the MIPS 70 fim im- 
age, the 850 yum SCUBA image shows an extended emission 
patch covering the entire area, in which a strong peak coincident 
with source A clearly stands out; 850 /im fluxes for the various 
24 jim sources were assigned following the same approach as 
for the 70 /urn fluxes. The peak at source A is also exactly coin- 
cident with the core visible in the interferometric OVRO images 
at 3.4 mm; the slight off'set of the 70;um peak position with re- 
spect to the 24yum source A-peak may then not be real and could 
be ascribed to an inaccuracy of the lOjjm image astrometry. 

Concerning the shorter wavelengths, for most of the sources 
we could estimate flux densities at 6.75 and 15 yum from 
ISOCAM images (Molinari et al. Il998b 1) using the same ap- 
proach followed for the MIPS 24 fim image; sources G, L, N and 
R are outside the ISOCAM FOV. Near-IR counterparts in J, H, 
and Kj could also be assigned for most of the sources (Faustini 
et al., A&A submitted). The K^image is shown in greyscale in 
the right panel of Fig. [T] the full-line contours report the 24 fim 
emission, while the dashed contours mark the emission pattern 
in the 3.6 cm thermal free-free emission (Molinari et al. |2002| l. 
The only dubious counterpart assignment concerns source H, for 
which there are three NIR objects very close to the 24 yum peak; 
we list all of them in Table [1] 

For those 24 fim sources with NIR counterpart we tried to ex- 
tend the SED coverage to the visible range by visual inspection 
of the Digital Sky Survey 2 blue and red plates. Where coun- 
terparts could be found, fluxes were estimated doing standard 
aperture photometry and calibrating the integrated flux using a 



conversion factor estimated from the magnitude of medium in- 
tensity sources in the field as reported in the USNO-Bl cata- 
logue. Zero-magnitude fluxes of 4260 Jy and 3080 Jy for the B 
and R bands were used (Bessel [T979l l. 

2.2. IRS Spectroscopy 

The areas observed in the 3x2 spectroscopic map are diff^er- 
ent across the wavelength range since the slits in the different 
modules have different sizes and orientations; however they are 
much smaller than the area imaged by MIPS. All slits were cen- 
tered on the millimeter peak position (coincident with source A), 
so that we have complete spectral coverage for this source. We 
used the CUBISM software (Smith et al. 2007) to mosaic and 
co-add the spectra from all the observed positions and obtain 
convenient spectral cubes for further analysis. We present in Fig. 
|2]the complete spectrum in surface brightness units extracted at 
the position of source A (top panel). 

Spectra from the two SLl and SL2 modules for A <15jum 
have been merged together and clearly show a set of Poly- 
Aromatic Hydrocarbon (PAH) features at wavelengths between 
6 and 13 yum. The SH spectrum is particularly noisy not only 
with respect to the lower resolution SL portion, but also with re- 
spect to the LH. Besides, it is also apparent (see bottom panel 
of Fig. |2]l that the SL and LH portions appear to be quite well 
aligned among them, while the SH portion seems anomalous in 
this respect. To understand the nature of the displacement of SH 
we resampled it to the same resolution as the SL portion, with 
which SH overlaps between 10 and 15yum; the resampled SH is 
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Fig. 3. The six identified lines in the entire spectram at the position of source A, as extracted from the spectrum in Fig.|2l the various 
lines are identified in each panel. 



Table 2. Integrated line fluxes at the position of source A. 



Line Wavelength Flux 



ipim) W cm - sr ' 



H2 (O-O)S(O) 


28.2 


6.4 X 10^ 


13 


H2 (O-O)S(l) 


17.02 


1.7 X 10- 


12 


H2 (0-0)S(2) 


12.28 


3.2 X 10" 


12 


[Nen] 


12.8 


3.2 X 10" 


12 


[Fen] 


26 


3.5 X 10- 


13 


[Sin] 


34.8 


4.3 X 10- 


12 



shown in black in the bottom panel of Fig. |2l normalising it to 
the SL over the entire overlapping range we obtain a ratio of 4.2, 
and if we rescale SH by this factor (the grey line in the figure) 
we see that SH almost perfectly overlaps with SL. We conclude 
that the SH spectrum likely suffers from an incorrect calibration, 
and a "gain" type of correction brings it in excellent agreement 
with both SL and LH. We then rescale SH by this gain factor, and 
we will use this rescaled spectrum to extract lines and estimate 
integrated fluxes. 

Figure [3] shows six lines identified in the SH and LH spec- 
tra. [Neil] 12.8//m line and the S(l) and S(2) pure rotational lines 
of H2 are extracted from the SH spectrum; the LH spectrum 
is relatively less noisy and allows us to reliably identify the 
[Fen]26yL(m, the H2 S(0) and the [Siii]34.8/im lines. Integrated 
fluxes were derived using Gaussian fits to line profiles after 
subtraction of linear baselines; estimated values are reported in 
Table 2 for the position of source A. 



3. Source A: an accreting pre-ZAMS massive YSO 

Source A is the most noticeable object in this region. While un- 
detected in the near-IR and ISOCAM images, it is first visible 
at 24 yum and then becomes the most prominent source in the 
far-IR, the submillimeter and the millimeter regions, and it co- 
incides with the location of the driving source of the molecular 
outfl ow seen in HC0+(1 -» 0) and SiO(v=0, 2 -> 1) (Molinari 
etaL ll998bl l. 

Source A, however, is not as prominent in the far-IR (Fig.|4]i 
as we originally believed when only IRAS data were available 
(Fontani et al. i2004l l. This has very important consequences re- 
garding its nature. 

In a systematic analysis of the SEDs of a sample of 42 
sites of massive star formation (Molinari et al. 2008 ) we used 
the radiative transfer models of Whitney et al. (2003) to com- 
pute a grid of SEDs for embedded ZAMS stars of various spec- 
tral classes (from B5 to 03) and for a wide range of enve- 
lope properties (radius, mass, geometry). In that work the far- 
IR portion of the SED was based on IRAS fluxes (the dia- 
monds in Fig. |4]i, corrected for a possible contribution from 
extended emission. We could not provide a successful fit with 
any of the models of our grid for what we now call source A 
in Moll60/IRAS23385+6053. Likewise, using the much more 
reliable estimate of the source SED in the far-IR range that the 
MIPS data allow (the squares in Fig. HI, we still cannot find any 
acceptable fit with models of embedded ZAMS stars; any model 
that fits the 850 /vm and 3.4 mm points predicts far too much flux 
shortward of 100 /im. We neglected the possible contribution of 
free-free emission to the 3.4mm flux as this could amount to a 
10-20% level at most, given that source A was not detected at 
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resentative of the region of the parameter space identified by the 
SED fitting. 

The greybody fit represented by the full line in Fig. |4] pro- 
vides an integrated bolometric luminosity of 3170 at a dis- 
tance of 4.9 kpc (similar to what was obtained by Stecklum et 
al. I2005I I. This is remarkably lower than the ~ 16 000 Lq ob- 
tained again fitting a greybody but using the IRAS PSC 60 and 
100 jum fluxes (the dashed line in Fig. [H see also Molinari et 
al. |1998bl ). or than the -10000 Lq obtained if we ti-y to correct 
the far-IR fluxes for extended emission using higher spatial res- 
olution information from the submillimeter, as done in Molinari 
et al. (2008). Thanks to the MIPS data we can now set much 
more stringent constraints on the SED over the entire wavelength 
range and we can obtain this much lower value for the luminosity 
while the estimated envelope mass remains basically the same. 



100 1000 
Wovelength (/xm) 



Fig. 4. Spectral Energy Distribution (SED) for the central core of 
Mo1160/IRAS23385-h6053. The diamonds are the IRAS fluxes, 
the squares are the MIPS fluxes and the asterisks represent the 
ensemble of measurements from ISOCAM (upper limits since 
there is no detection below 24 /jm), SCUBA and OVRO; the 
dashed line is a greybody fit to the above data set using the IRAS 
60 and 100 //mfluxes (as in Molinari et al. 1998b). The full Une 
is a greybody fit to the SED using the MIPS data instead of the 
IRAS ones. 



2cm at the VLA in the B-configuration (Molinari et al. |1998a] l. 
and in the extreme case of optically thick free-free (Foe v^) up to 
the milhmeter range. 

The full line in Fig. |4] represents the best fit, using a simple 
greybody which assumes constant density and temperature. The 
fit uses the OVRO, SCUBA, and MIPS data, and the ISOCAM 
upper limits (the asterisks) but neglects the 12 and 25 fim IRAS 
data since they are clearly due to the other YSOs in the region 
(see Fig.[T]). The overall fit is remarkably good and corresponds 
to an envelope with T=37 K, M=220 Mq and a dust opacity 
index jS=l, assuming a dust opacity of 0.005 cm^g"' at 1.2 mm 
and a gas-to-dust ratio of 100 (Preibisch et al. |1993l ). 

Robitaille et al. (2006) used the same radiative transfer code 
of Whitney et al. (2003) that we use in Molinari et al. (20081), 
to produce a much more extensive grid of SED models (not yet 
available at the time of our previous work) including central stars 
in all pre-MS phases from the birthline to the ZAMS. It is quite 
convenient to use the automatic SED fitting tool provided by 
Robitaille et al. ( 120071 1, which reports the best fitting models in 
increasing;jf^ order. Fitting the SED of source A with this grid of 
models provides a best fit for a central object of ~9.5 M© with a 
radius of ~25 R© and a surface temperature slightly in excess of 
8000 K accreting at ~ 2-10"^ Mq yr"' from a massive ~ 500 Mq 
envelope. This value of the stellar radius is much higher than 
typical values for ZAMS stars of that mass; besides, with these 
accretion rates the stellar mass needed to ignite deuterium burn- 
ing is higher than 20 Mq (Palla & Stahler 1992), so that such 
an object would still be in the pre-ZAMS phase. Clearly, given 
the number of parameters involved in the modelling and the lim- 
ited number of data points available, the fitting cannot pinpoint a 
unique combination of parameters but rather identifies a certain 
range in the parameter space. It is true, however, that the possi- 
bility of an embedded ZAMS star seems to be reliably excluded 
and that the above-mentioned combination of parameters is rep- 



10* 

10^ 
10^ 
10' 
10° 




100 



Fig. 5. Lboi-Menv diagram from Molinari et al. ( 120081) : full circles 
are IR-P sources, empty circles are MM-P sources and plusses 
are for IR-S sources. Dotted full lines and diamond dashed lines 
are evolutionary models; the two squares connected by the ver- 
tical dotted line represent the extent of the shift of source A due 
to the new luminosity estimate. The full and dashed diagonal 
lines represent a linear fit to the position of the IR-P and MM-P 
sources, respectively. 

This decrease in the luminosity estimate has important con- 
sequences concerning the evolutionary stage of the massive 
forming object. In our recent analysis of the SED evolution of 
massive YSOs (Molinari et al. 2008 ) we could differentiate two 
classes of objects depending on whether the entire SED from the 
mid-IR to the millimeter can be consistently fitted with a single 
model of an embedded ZAMS star (called "IR-P" by Molinari 
et al. and represented with the filled circles in Fig. |5]) or not. 
If not, the SED was modeled with two components: a single- 
temperature greybody for the millimetric/far-IR part ("MM-P", 
open circles in Fig.|5]l, and a mildly-obscured ZAMS star for the 
mid-IR ("IR-S", plus signs in Fig.|5]l. These different classes of 
objects occupy distinct regions in the Lboi-Menv diagram (Fig.|5]l, 
and a simple toy-model for the evolution of the bolometric lumi- 
nosity and the circumstellar envelope can explain the sequence 
MM-P/IR-P/IR-S in evolutive terms analogous to the Class O-I- 
II sequence established for low-mass YSOs. From a statistical 
viewpoint the IR-P objects are dominated by sources with far-IR 
colours of Hot Cores/UCHii regions, therefore confirming their 
ZAMS nature; this is not the case for MM-P objects. In addi- 



Molinari et al.: Pre-ZAMS nature of Moll60/IRAS23385+6053 



7 



tion to the SED differences, Fig.|5]shows that the luminosity is a 
critical parameter that can trace with a high dynamical range the 
transition between MM-P and IR-P objects. 

The two squares in Fig.|5]illustrate the position of the source 
with the luminosity that was found in our original work on this 
source (Molinari et al. |1998bT l and the one determined using the 
MIPS observations. Quite remarkably, the shift to lower lumi- 
nosity brings the source from the locus of the IR-P objects (the 
full straight line, which in our interpretation marks the arrival 
onto the ZAMS) to the locus of the MM-P objects (the dashed 
straight line). In evolutionary terms within the framework of the 
model proposed in Molinari et al. (120081) . this luminosity cor- 
rection makes the source more than 5 10"^ years younger than 
we have previously estimated; this is the time it takes, according 
to the prescriptions of accelerating accretion models (McKee & 
Tan 2003) to bring source A onto the ZAMS from its present 
location in the Lboi-Menv diagram. Source A is even younger if 
we adopt the ~ 500 Mq envelope mass derived from the detailed 
grid of SED models of Robitaille et al. ( |2006l ). 

Molinari et al. (2008) also suggested that MM-P objects are 
characterized by a much steeper SED for A < IQO fim compared 
to IR-P objects. In particular we predicted that the [24-70] colour 
(estimated by interpolating observed mid-IR and sub-mm spec- 
tra using SED modeling) to be ~1 for IR-P objects and ~4 for 
MM-P objects: the MIPS observations provide a direct measure- 
ment of the SED shape in this critical region: with [24-70] =2.6 
for source A, these measurements strenghten the foundations for 
the working assumption in our previous work. 

4. Other YSOs and the UV-field in the region 

Since source A is only detected longward of 24 fim in the con- 
tinuum, dust extinction will prevent any line emission at shorter 
wavelengths, originating from the embedded YSO, to emerge 
and be detectable. Therefore all PAH features, the [Neii] 12.8/vm, 
and the two S(2) and S(l) lines cannot originate from the mas- 
sive central core; indeed all emission features are detected with 
variable intensities throughout the regions mapped with the IRS. 
To better investigate their spatial distribution we again used 
CUBISM to produce continuum-subtracted spatial maps of the 
integrated line fluxes. They are presented in Fig. |6] for the 
11.3 yum PAH feature and for five of the detected lines, as in- 
dicated. 

The intensity of the 11. 3 fim feature (the same spatial distri- 
bution is followed by all other PAHs) has a minimum at the lo- 
cation of source A and increases towards the north where source 
F is located. PAHs emit in a non-equilibrium regime in the pres- 
ence of relatively energetic incident radiation field. Their detec- 
tion, with integrated intensities similar to that of the underly- 
ing continuum, suggests then the presence of an irradiating UV- 
field in the region. Indeed, all detected transitions are commonly 
found in regions of intense UV-irradiation where the local ISM 
is photo-dissociated ([Fen]26/im, [Siii]34.8;um plus the three H2 
lines; see Tielens & Hollenbach 1985 and Kaufman et al. 1 19991 1 
or photo-ionized ([Neii]12.8/im, see Morisset et al. 120021 ). con- 
sistent with the spatial distribution of the continuum-subtracted 
integrated emission for all detected lines. In all panels of Fig.|6] 
the grey-scale shows, with the exception of H2 S(0), that the 
intensity of the continuum-subtracted integrated line fluxes in- 
creases away from source A. We will discuss this in more detail 
below, but we anticipate that this distribution seems to rule out 
any significant role of source A in the budget of the UV-field in 
this region. The H2 rotational lines are more commonly inter- 
preted as being due to C-shocks, but the extended distribution 



observed seems more consistent with a PDR origin. The line ra- 
tio analysis using all the detected lines in Fig. [3] was carried out 
using the PDR Toolbo?|3 that conveniently implements a web- 
based fitting program using the diagnostic diagrams based on 
the recent models by Kaufman, Wolfire & Hollenbach (I1999I I. 
Metallicities 2-3 times higher than local ISM values, and more 
typical of Hii regions and massive star forming regions in gen- 
eral (Shaver et al. 1983; Tielens & Hollenbach 1985 ) have been 
used. Most of the line ratios are consistent with a far-UV field 
intensity of the order of 100< Go < 1000, in units of the FUV 
field intensity in the solar neighbourhood (Habing 119681 ), for 
an ISM density in excess of 10"* cm""*. The lines which deviate 
most from defining a consistent picture are the H2 S(2) and the 
[Feii]26//m, which appear to be a factor of 2-3 stronger than PDR 
models predictions given the fluxes of the other lines. The mod- 
eled ISM densities are entirely reasonable for the intra-clump 
material of a massive star forming region, as shown from maps 
in high-density tracing species like HCO^(l-O), '^CO(2-l) and 
CS(3-2) (Brand et al. l20UTT ). 

To verify if external field objects could possibly be respon- 
sible for the detected UV-irradiation, we searched the Spectral 
Classification catalogues available in Aladin for stars with spec- 
tral type earlier than B. The nearest object is the B9V star 
BD-H60 2600 located more than 10' to the North; assuming stan- 
dard stellar parameters (Thompson i 1 984i l and a distance equal 
to Moll60/IRAS23385+6053 the resulting Go is few orders of 
magnitudes below the observed levels. 

Among the various detected sources in this region, source A 
is the most massive and luminous; its extreme optical depth be- 
low 20jum, however, would in principle allow only [Feii]26yum, 
[Siii]34.8jum and H2 S(0) to escape from its inner envelope re- 
gions. However, the only indication we have in this respect, as 
already noted, is limited to the H2 S(0) which however traces 
warm gas in general; the former two lines, that are the most 
commonly used PDR tracers, are instead peaking away from 
source A. The fact that also the intensity of the PAH feature, 
which is directly related to the intensity of the UV-field, in- 
creases away from source A confirms that the latter cannot be 
the radiating source of the UV-field, not even in the presence of 
a clumpy circumstellar envelope. 

The other sources detected in the region are much brighter 
than A at 24 /urn and many have a near-IR counterpart, so that 
they could be intermediate mass YSOs in a more evolved state 
than source A. However, given the limited spatial coverage of 
the IRS maps, it is not possible to use spectroscopic tools to get 
more insight into the nature of these various objects. 

More can be learned from the SED modelling of the other 
sources revealed at 24//m (Table [1). Now, we seek evidence 
for the presence of an embedded ZAMS star, as opposed to a 
pre-MS object, which might be responsible for the relatively in- 
tense UV-irradiation conditions in the region, as well as for the 
extended radio continuum which is thermal free-free in origin 
(see Fig. [T}?). We again used the automatic SED fitting tool of 
Robitaille et al. ( 120071 ) . For sources with no counterpart in the 
near-IR or in the visible, we adopted the limiting magnitudes 
of our own JHKj data (Faustini et al., submitted) and of the 
DSS2 plates. We searched, among the various ^^-ranked fits, 
for models where the central stars have photospheric tempera- 
tures and stellar radii typical of ZAMS stars, to verify if some 
of the sources detected in the region could be capable of radi- 
ating sufficient quantities of UV flux to photo-dissociate the in- 
traclump medium as suggested by the IRS spectroscopy. Such 



' The PDR Toolbox is at http://dustem.astro.umd.edu/pdrt/index.html] 



Molinari et al.: Pre-ZAMS nature of Moll60/IRAS23385+6053 



61°10'40" 



61°10'20" 



o 
o 

^61°10'00' 

"-a 



61°10'40" 



61°10'20" 



61°10'00" 




23''40'"55^ 



23*'40™55^ 



50^ 23''40'"55^ 



a(J2000) 



Fig. 6. Grey-scale maps (intensity increases from white to black) of continuum-subtracted integrated line intensities superimposed 
on the the contours of the 24 //m continuum emission. The various panels show, from left to right and from top to bottom: the PAH 
11.3 jum feature, H2 S(2), H2 S(l), [Feii]26//m, H2 S(0), and [Siii]34.8//m. The map for the [Neii]12.8yum line is not reported but the 
emission distribution is vutually identical to the H2 S(l) and S(2) lines. 



models could be found for sources B, C, D, H and I, although 
not with the best formal they rank lower in the best list, 
but still provide very good fits (shown in Fig. |7]l. No such fits 
could instead be found for the other sources in Table [1] which 
can only be fitted with pre-ZAMS objects. Table [3] reports the 
model parameters for SED of embedded ZAMS stars; Col. 9 
lists the Lyman continuum as tabulated by Thompson (,1984i) for 
ZAMS stars with photospheric temperatures similar to those re- 
ported in Col. 3. Col. 8 reports the model intrinsic bolometric 
luminosity. The sum of the values in Col. 8 plus the luminosity 
previously estimated for source A amounts to ~30 OOOL0 which 
is 50% higher than the value obtained for the entire region using 
the IRAS fluxes; we regard this discrepancy as barely significant, 
given the spread of the model luminosity depending on the exact 
model selected in the fit. Col. 10 shows the unattenuated inte- 
grated 6-13.6 eV far-UV continuum, expressed in units of the 
far-UV field intensity in the solar neighborhood (Habing[T968|l, 
irradiated by each source at the position of source A. 

The total Gq at the position of source A is consistent with the 
regime deduced from the infrared IRS lines, if PDR in origin. 
Since Go in Col. 10 of Table [3] is computed with no interven- 
ing attenuation between the emitting source and the position of 
source A, we must verify that this is indeed plausible. Indeed, all 
the models in Tableware obtained for systems where the ZAMS 
star is surrounded by a disk with high inclination angle with re- 
spect to the line-of-sight, and with no envelope. The SEDs which 
fit the observed data points (see Fig. |7]i are obtained assuming 
an additional extinction correction along the line-of-sight to the 
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observer (Col. 7) due to intervening dust not related to the im- 
mediate circumstellar environments. Since these values are not 
homogeneous, and the distribution of optically visible objects in 
the Moll60/IRAS23385-i-6053 area is not suggestive of signifi- 
cant or variable interstellar absorption, we believe this extinction 
to be intra-cluster in origin. It is clear that with such amounts of 
extinction the far-UV field should be entirely absorbed, so that 
we have to assume that the line-of-sight from the various sources 
to source A must be relatively dust-free. That this must be the 
case also seems to be suggested by the distribution of the ra- 
dio continuum from thermal free-free emission (dashed lines in 
Fig- [It); the extent of the emission requires that sources of ion- 
izing continuum must be present, and that this continuum is not 
confined to the immediate surroundings of the ionizing stars. The 
most likely sources for the radio emission seem to be sources D 
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Fig. 7. SEDs for the various sources where a fit with a model of an embedded ZAMS star could be found. 



and H for the eastern radio lobe, and sources B and C for the 
western radio lobe. 

Concerning the intensity of the radio emission, however, its 
conversion into unattenuated Lyman continuum intensity N^y 
provides values (Molinari et al. 120021 ) one order of magnitude 
higher than those that can be justified by the modeled ZAMS 
objects in the region (Col. 9 of Table [3]). Even integrating the 
IRAS and millimeter fluxes would yield a value of ~20 000 L© 
which, also if coming from a single ZAMS star (which is not 
the present case), is again not sufficient to provide the amount 
of Niv deduced from the radio fluxes. It should also be noted 
that the latter is a lower limit because the free-free emission has 
been assumed optically thin and possible attenuation from dust 
inside the HII regions has not been considered. The simplest 
way to reconcile this apparent inconsistency would be to posit 
that the distance of 4.9 kpc which we assumed for this object 
is incorrect. Indeed, while the bolometric luminosity scales with 
the square of the distance, the relationship between ZAMS lu- 
minosity and Lyman continuum is much steeper In the limiting 
case of a single object emitting all the luminosity (which is un- 
realistic given the evidence), moving the object to a distance of 
~6 kpc would raise the luminosity to ~29 000 L© which would 
produce Log(N£v)~47.4; this is equal to the Nl^ deduced from 
the observed radio flux rescaled to a distance of 6 kpc. The ev- 
idence from the radio continuum, however, is that there are at 
least two major sources of the ionizing-flux field; a simple cal- 
culation shows that in the limit of two dominant ZAMS objects 
the distance should be increased to about 8 kpc. The "Vlsr that 
would be implied by this distance is more than 20 km s ' differ- 
ent from the value measured from radio spectroscopy, and this is 
not incompatible with the magnitude of the streaming motions 
deduced from the observed velocity field (Brand & Blitz |1993l ). 



Changing the distance of this region to 8 kpc does not af- 
fect our interpretation about the nature of source A as it emerges 
from Fig. |5j both envelope mass and luminosity scale the same 
way with distance so that the position of source A in Fig. |5] 
would shift toward upper-right along a line of slope 1, there- 
fore leaving the source in the area occupied by MM-P objects. 
This higher distance would also imply that the mass regime of 
the Moll60/IRAS23385+6053 region is higher than previously 
thought; in particular source A, which in our proposed evolution- 
ary interpretation and at the distance of 4.9 kpc would reach the 
ZAMS (the solid line corresponding to the location of the IR-P 
objects in Fig. |5]) at a luminosity of nearly I6OOOL0 as a BO. 5 
star, if placed at the distance of 8 kpc would reach the ZAMS at 
about 40 OOOLq, roughly coiTesponding to an 09.5 star. 



5. Conclusions 

Sub-arcminute spatial resolution in the far-IR proves to be crit- 
ical for a proper assessment of the evolutionary stage of a mas- 
sive forming object. SED modeling of the most prominent YSOs 
in the Moll60/IRAS23385+6053 region suggests the presence 
of several intermediate and high-mass YSOs in different evo- 
lutionary stages. While some of them are compatible with be- 
ing ZAMS objects with spectral types between BL5 and B5, 
source A is best interpreted as a strongly accreting object of 
central mass comparable to the other ZAMS YSOs in the area, 
but not yet on the ZAMS. The region immediately surrounding 
source A is bright in infrared lines which are typical of photo- 
dissociation/ionisation regions; the deduced intensity of the ra- 
diative far-UV field at the position of source A is compatible 
with emission from the few ZAMS objects as characterized by 
the SED modelling, which are also likely to be responsible for 
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the extended radio continuum emission. There is an inconsis- 
tency, however, between the Lyman continuum estimated from 
the radio flux and from the luminosity of the modeled ZAMS 
stars in the region. At the moment we can only reconcile this if 
we assume that the region is almost twice as distant as assumed 
up to now, which obviously would need independent confirma- 
tion. 

The extreme values of the circumstellar extinction implied 
by the extreme SED of source A exclude that this source can 
participate in any way to the far-UV radiative field which per- 
meates the region. As a whole, the population of intermediate 
and massive YSOs in the region is suggestive of a star formation 
timescale of the order of few 10^ years, as suggested by models. 

However, intermediate- and high-mass stars are not the only 
objects found in this region. A cluster of lower-mass stars is re- 
vealed in the near-IR; this is partially apparent in Fig.[TJ), where 
the contrast, however, is not optimal to visually reveal the much 
more abundant population of fainter objects. Stellar density anal- 
ysis shows the cluster to be as extended as the extent of the mid- 
and far-IR emission, with several tens of members (Faustini et 
al., submitted). The comparison of the K-band luminosity func- 
tion with an extensive grid of synthetic cluster models obtained 
for a wide range of stellar ages, IMFs and star formation his- 
tories, suggests that the ages of these lower mass objects are at 
least few 10* years. 

The simultaneous presence of relatively old pre-main se- 
quence objects with a massive YSO still in an active, pre-ZAMS, 
accretion phase, confirms that star formation tends to be an on- 
going process for quite a long time span. The highest-mass star 
seems to be the last one to form. 

The numerical importance of massive YSOs in a pre-ZAMS 
phase cannot be firmly established at the moment. Submillimeter 
surveys (Hill et al. 2005 , Beltran et al. 12006^ reveal a consis- 
tent population of submillimeter cores in the proximity of IRAS 
point sources. These cores are devoid of mid-lR emission, but 
without an assessment of their far-infrared properties which only 
will allow firm temperature, luminosity and mass estimates, their 
nature will remain elusive. The PACS and SPIRE far-infrared 
cameras on board the Herschel satellite will be the ultimate 
tools to help distinguish between Hot-Cores, pre-ZAMS objects 
and quiescent cold cores. The recently approved Herschel Key- 
Project Hi-GAL for a complete 60-600/im continuum survey of 
the inner Galactic Plane will be invaluable to obtain firm statis- 
tics of massive YSOs in all evolutionary stages, thus providing 
a solid foundation to establish the timeline for the formation of 
intermediate- and high-mass stars. 
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